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Abstract. Remote observations of coronal holes have strongly implicated the kinetic interactions of 
ion-cyclotron waves with ions as a principal mechanism for plasma heating and acceleration of the 
fast solar wind. In order to study these waves, a linear perturbation analysis is used in the work frame 
of the collisionless multi-fluid model. We consider a non-uniform background plasma describing a 
funnel region and use the ray tracing equations to compute the ray path of the waves as well as the 
spatial variation of their properties. 



1 Introduction 



The ultraviolet spectroscopic observations made by SUMER and UVCS aboard SOHO in- 
dicated th at heavy ions in the coronal holes are very hot with high temperature anisotropy 
(see, e.g. jKohl et al.lll99"7Hwilh elm et al. 1998). This result is a strong indication for heat- 



ing by ion-cyclo tron resona nce (i.e., collisionless energy exchange between ions and wave 
fluctuations, see IStixiri99 2'. Chap 10) involving ion-cyclotron waves that are presumably 



generated in the lower corona from small-scale reconnection events ( Axford & McKenzj^ 



1995|) . Performing a Fourier plane wave analysis, ion-cyclotron waves are studied using the 
collisionless multi-fluid model. While neglecting the electron inertia, this model permits the 
consideration of ion-cyclotron wave efl'ects that are absent from the one-fluid MHD model. 
Realistic models of density and temperature as well as a 2D funnel model describing the 
open magnetic field are used to define the background plasma. Considering the WKB ap- 
proximation, we first solve locally the dispersion relation and then perform a non-local wave 
analysis using the ray-tracing theory, which allows to compute the ray path of the waves in 
the funnel as well as the spatial variation of their properties. 

2 Basic equations 

The cold collisionless fluid equations for a particle species j are: 

-l + ^.(„.y.) = 0, (1) 
""PA-^ + ■ "^^i) - W(E + vy X B) = , (2) 
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Figure 1. Left: El ectronic density (Ue) and temperature (T) model-profiles of the chromosphere 
Jpontenla et alJll993h and the lower corona (Gabriel 1 97^ . Ri^ht: Funnel magnetic field geometry 
as obtained from a 2-D potential field model (IHackenberg et alJuOOd) . 



where mj, rij and are respectively the mass, density and velocity of a species j. The 
electric field E and the magnetic field B are given by Earaday's law, i.e. V x E = -dB/dt. 
The background density, temperature and magnetic field are described in Eig.[T] Considering 
the quasi-neutrality and no ambient electric field, we perform a Eourier plane- wave analysis. 
The dispersion relation is obtained using 

D(oj, k, r) = Det[—k x (k x E) -h s(aj, k, r) • E] = , (3) 

where c is the speed of light in vacuum, s is the dielectric tensor which is a function of the 
wave frequency oj, the wave vector k and the large-scale position vector r. The wave vector 
k lies in the x - z plane, with k = k(sin6, 0, cosO). 

The dispersion relation © is used to compute the Hamiltonian-type ray-tracing equations 
dWeinbergll 19621) : 

dr _ dDjoj, k, r)/dk dk _ dDjco, k, r)/dr 

dt ~ ~dD(oj, k, r)/doj' 'dt ~ dD(oj, k, r)/doj ' ^ ^ 

3 Numerical results 

We consider a cold plasma made of electrons (n^), protons (rip) and alpha particles (He^^, 
indicated by a) with ria = O.lrip. In this case, the dispersion relation 0]) is a quadratic 
polynomial of degree 6, which means that 3 modes exist and each one is represented by 
an oppositely propagating (oj > and oj < 0) pair of waves. These modes, for which the 
dispersion curves are shown in the right panel of Eig. [21 are the ion-cyclotron modes 1 (ICl) 
and 2 (IC2) and the fast mode. Eor large k, the two IC modes reach a resonance regime at 
each one of the ion-cyclotron frequencies, i.e. oj = flp and oj = fla = The fast mode 

has a cut-ofif frequency ojcol^p = 0.583 and couples and mode converts with the ICl mode at 
the so-called cross-over frequency ojcrl^p = 0.612. The polar plots of the group velocity for 
kVAp/^p=0.2 (Eig.O, show that the IC2 mode is mainly anisotropic and cannot propagate 
perpendicularly to Bq. Thus the corresponding energy mainly flows along the magnetic field 
lines, with maximum group velocity at parallel propagation. On the other hand, the ICl 
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Figure 2. Right: Dispersion curves in the cold three-fluid (e-p-He^+) model at the location (x=7.5 
Mm, z=2.5 Mm) in the funnel with a Bq -inclination angle (p ^ 79° and for different propagation angles 
6 (Qp and V^p are the proton cyclotron frequency and Alfven speed). Left: Polar plot of the group 
velocity as a function of 6 and for kVApl^p - 0.2. 



mode and the fast mode propagate almost isotropically, and consequently the energy flow is 
fairly isotropic. These results are confirmed by the ray-tracing computation which clearly 
shows that the IC2 mode is well guided along the field lines, since the ray paths (direction 
of the group velocity) for various initial angles of propagation nicely follow the magnetic 
field lines (Fig. [3]). Indeed, the maximal angular deviation is</^^6.5°atz^^4 Mm and 
i// ^ 0° (quasi-parallel propagation) in the upper part of the funnel, where the wave is quasi- 
electrostatic, i.e. I ^ |~ 1, and has a nearly linear polarization, i.e. | ^ |~ 0. On the contrary to 
that behavior, the ICl mode and the fast mode are unguided with a ray path having mainly a 
strait trajectory. The angular deviation i// between the ray path and Bq varies between -60° 
and 80° for the ICl mode and -40° and 60° for the fast mode. The polarization is in general 
elliptical, right-handed (g > 0) for the ICl mode and left-handed (g < 0) for the fast mode, 
except for = 20° in the upper part of the funnel (z> 6 Mm) where it is mainly circular, 
right-handed (^ = -1) for the ICl mode and left-handed (^ = -1) for the fast mode. 
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Figure 3. Ray paths of the waves and the spatial variation of their basic properties. The waves are 
launched at the location (xo=7.5 Mm, zo=2.5 Mm) in the coronal funnel (with a Bq -inclination angle 
^0 ~ 79°) with an initial normalized wave number ko=0.2 and different initial angles of propagation 
^0- In the top panels, the dashed lines represent the funnel field lines. The wave properties are: {g) 
helicity (degree of circular polarization), {if/) angle between the direction of the group velocity and Bq, 
(^) the electrostatic part of the wave. 



